o 
o 
o 

(N 
•*-> 
o 

O 
o 


> 
o 
l> 
o 
o 

o 
o 


Oh 


X 


Comment on "First Order Transition in the 
Ginzburg-Landau Model" 

In a recent Letter, Curty and Beck jjj have shown 
very interesting results which indicate that the Ginzburg- 
Landau (G-L) transition becomes first order when the 
coherence length £=£ |i|~2 (t = T/Tq-1 is the reduced 
temperature) becomes of the order of the lattice spacing 
e. They considered the lattice G-L model parametrized 

2 

by: a = which controlls the strength of amplitude 

fluctuations (that grow as er decreases), and Vo = j^-§7> 
which governs the overall strength of the complex G-L field 
tj) = \~ip\exp[i9]. They treated the model by a variational 
approximation and got a criterion for first order transition 
in the form of an inequality which involves the spatial di- 
mension d. They concluded that their criterion for d = 2 is 
not clear and that doubt remained if the first order found 
for this case is not an artifact of the used approximation. 

Here I will present clear evidences that for d=2 a first 
order transition takes place when £ becomes ~ e {a ~ 1) 
and that this is connected with a sudden proliferation of 
vortices. Similar results where reported in Q although 
using a different parametrization of the G-L model which 
obscures the comparison with 0. The d — 2 G-L hamil- 
tonian H [a, Vo) was simulated on L x L lattices. The mea- 
sured phase diagram on the plane T-d~? is showed in 
Fig.I where, as in ref. §J|, f = T/(-tV ) and a = -to. 



Figure I: Phase diagram in the (T, a 2 ) plane for L = 40. 

Above £/a = a~i = 0.8 the phase transition line 
changes from second order (dashed line) to first order 
(filled line). The double peak of the energy density e his- 
togram corresponding to the two coexisting phases, char- 
acteristic of a first-order transition, is showed in Fig. 2- 
a for <t - 2 = 0.85. Both peaks remain fixed as L in- 
creases and the width of each of them clearly scales as 
x/Gjtt) = due to ordinary non-critical fluctuations 
(Fig. 2-b). In addition, a strong hysteresis effect was 
found for e when considering heating and cooling runs. 
On the other hand, for <r - 2 < 0.75 the peaks are much 
lower and wider, they move towards to an intermediate 
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Figure 2: Histograms of e for L=10 (below), L—20 (mid- 
dle) and i=40 (above), a) <r~2 = 0.85. b) Zoom of the 
right peak, c) <7~2 — 0.75. 

value of e as L increases and their width do not scale as 
■£ (Fig. 2-c). Moreover, no hysteresis in e is found. 

The central role played by vortex excitations in de- 
termining the nature of the phase transition can be seen 
in Fig. 3-a where v is plotted vs. T for different values 
of a and L=40. For a—1 there is a clear discontinuity in 
the vortex density v (Fig. 3-b) . As long as a increases the 
jump becomes more smooth and moves to higher values 
of T c until for a = 100 one gets something very close to 
the Kosterlitz Thouless (K-T) behavior of the XY model. 
The enhancement of vortex production when amplitude 
fluctuations are large is due basically to the fact that they 
decrease the energy of vortices. 



Figure 3: a) v vs. T for a=l (A), or=1.5625 (V), ct=4 
(x), ct=100 (+), XY model(o). b) Zoom of g-(a) 

Therefore, in the G-L model the nature of the phase 
transition depends dramatically on the value of a: For 
a <1.5625 (- > 0.8) the density of vortices experiments 
a discontinuous jump which coincides with a first order 
transition. On the other hand, for a ^> 1 the G-L re- 
duces to the XY model with the more subtle K-T phase 
transition. 
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